Magnetization in Molecular Iron Rings 
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The organometallic ring molecules Fe6 and Feio are leading examples of a class of nanoscopic molec- 
ular magnets, which have been of intense recent interest both for their intrinsic magnetic properties 
and as candidates for the observation of macroscopic quantum coherent phenomena. Torque mag- 
netometry experiments have been performed to measure the magnetization in single crystals of both 
systems. We provide a detailed interpretation of these results, with a view to full characterization 
of the material parameters. We present both the most accurate numerical simulations performed 
to date for ring molecules, using Exact Diagonalization and Density Matrix Renormalization Group 
techniques, and a semiclassical description for purposes of comparison. The results permit quantita- 
tive analysis of the variation of critical fields with angle, of the nature and height of magnetization 
and torque steps, and of the width and rounding of the plateau regions in both quantities. 

PACS numbers: 75.10.Jm, 03.65. Sq, 73.40. Gk, 75.30. Gw 
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I. INTRODUCTION 

The molecular iron ring, or "ferric wheel" systems Fe6 
and Feiocl present an interesting subgroup of magnetic 
molecular clusters. Both materials have dominant anti- 
ferromagnetic (AF) coupling between spins S = 5/2 on 
each iron site, and have a ground state of total spin S = 0. 
The magnetic properties of such nanoscopic molecules 
result from the interplay of superexchange interactions 
between the atomic spins, dipolar coupling of the local 
moments, and on-site spin anisotropies arising from lig- 
and configurations. In view of these complexities, and 
of the possibilities which existEl to influence the relative 
strengths of each interaction type through the non-Fe 
constituents of the molecules, a detailed understanding 
of the magnetic response is crucial. In appropriate pa- 
rameter regimes, the ring systems are consideredou to be 
candidates for the observation of "macroscopic" quantum 
phenomena, in the form of quantum coherent tunneling 
of the Neel vector. 

The magnetization M(B) of single crystals of both Fcg 
and Feio has recently been studied experimentally by 
cantilever torque magnetometryja using magnetic fields 
up to 23T, and at temperatures down to 0.45K. In the Fe6 
system, a selection of materials is available:El for Na:Fe@ 
only the first magnetization step, to the 5 = 1 state and 
occurring at the field B c \ ~ 17T, is accessible; for Li:Fe6, 
B C 2 ~ 21T is also observed. For Feio, the fields B c \ 
to B c 4, corresponding to the lowest four magnetization 
steps, may be measured. In a very recent analysis^ of 
Na:Fe6, temperatures as low as 30mK were probed, and 
a temperature-independence of the step width observed, 



a result to which we will return later. In each mate- 
rial, B c i is found to be largest in the orientation where 
the magnetic field is applied along the axes of the rings 
(which we will denote as ip = 0), and smallest when the 
field is directed in the plane of the rings (ip = tt/2). The 
variation between these extrema is to a good approxima- 
tion a sinusoidal function of 2ip. 

Here we present a detailed analysis of the magnetiza- 
tion in these ferric wheels, considering the spacing and 
widths of the magnetization steps, the flatness of the 
plateaus, and the angular variation of the critical fields. 
The primary results are given by numerical simulations: 
we perform both the first Lanczos diagonalization of the 
full Fe6 system, and the first Density Matrix Renormal- 
ization Group (DMRG) studies of finite, S — 5/2 sys- 
tems. Both techniques are exact at zero temperature 
for 6-site rings, and DMRG gives highly accurate results 
also for the 10-site case. As an aid to interpretation, 
we appeal in Sec. II to a semiclassical formulation for 
the molecular spin as a rigid-rotor model, and for the 
low-energy dynamics of the Neel vector.cl In Sec. Ill we 
compare this model with existing magnetization data to 
extract the relevant materials parameters. In Sec. IV we 
present the numerical simulations, and discuss the qual- 
itative and quantitative significance of the results. Sec. 
V contains a summary and conclusions. 



II. ANALYTICAL DESCRIPTION 

We begin by considering a semiclassical description for 
the behavior of magnetic molecular rings, in order to de- 
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duce approximate materials parameters for the simula- 
tions. Justification for the applicability of the semiclassi- 
cal model may be found in Ref. ||. The minimal starting 
Hamiltonian for the ring system is 
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h = jj2 Si-Si+i + ^(s*) + hh - E s » 



(1) 



1=1 



where N — 6 or 10, Si = Sjv+i, and h = j^ibB/S. J is 
the AF superexchange interaction, U contains intra-ring 
dipolar interactions and single-ion anisotropy terms, and 
the final term is the Zeeman coupling. For a ring geom- 
etry, it is easy to show that dipolar interactions within 
a single molecule contribute only an effective easy-axis 
spin anisotropy which favors the axis normal to the ring 
plane. In addition, the environment of each Fe ion due to 
the surrounding ligand groups in general has a non-cubic 
point-group symmetry; for the majority of cases the on- 
site symmetry may be taken as uniaxial, and, again with 
the assistance of the ring geometry which allowsu radial 
anisotropy to be subsumed within the uniaxial terms, we 
will assume this form of anisotropy. In all systems known 
to date, this has been found to be weak, and of easy-axis 
orientation in the same sense as above. Thus to a rea- 
sonable approximation, the separate contributions to the 
U term may be combined as a single, effective on-site 
anisotropy term X^i=i ~^zS 2 z , in which the sign of k z 
is chosen such that a positive value implies an easy axis 
for the orientation of the local spins. In Ref. || it was 
shown that the dipolar and single-site components of the 
effective k z may be attributed separately. 

The Hamiltonian (§) is mappedu to a form of nonlin- 
ear a model for the staggered magnetization n, by ex- 
pressing the spins Si in coherent state representation as 
Sfti = (—lySrii + li, and integrating over the uniform 
component 1. One obtains 

C E =N^[h 2 + »h.(n Ah) - (h 2 - (h • n) 2 )] 

~k z S 2 n 2 z ) . (2) 

Gradient terms (c^n) 2 , where x is the spatial coordinate 
around the ring, may be neglected for small N, because 
for a small, closed system such as the ring geometry they 
represent significantly higher- lying energy states. This 
semiclassical description of the low-energy physics,is also 
suitable for analyzing quantum coherence effects.u 

The Hamiltonian formulation canonically conjugate to 
Ce @j is deduced by considering the continuous spatial 
variable x = Tin, its conjugate momentum p = dCs/dic, 
and the angular momentum L = x A p. In real time 
t = —it this is 
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Nk z S 2 n 2 



(3) 



which may be split in obvious notation into Hq(L) + Ha, 
the latter term denoting the anisotropy contribution. Hq 



has eigenstates \L, L z ) — \l, m), with I = 0, 1, 2, . . . and 
— I < m < I, and eigenvalues Ei <m — (2J/N)l(l + 1) + 
guBBm. Thus Hq describes a staircase-like magneti- 
zation with steps occurring at the field values B cn — 
AJn/gnBN , where n = 0, 1, 2, ... is an integer. This be- 
havior accounts directly for the gross features of the ob- 
served magnetization curves.0 Inclusion of the anisotropy 
term Ha allows one to proceed to the effective spin 
HamiltonianLl employed in Ref. |^, and to a discussion 
of the variation of the critical field with angle, B cn (ip). 
In practice, this is used to deduce the anisotropy param- 
eter k z . In the presence of a single-ion anisotropy, the 
energy levels at zero field of the first (triplet) excitations 
are not degenerate, and B cl is given by the field where 
the m = — 1 level crosses the singlet. For arbitrary ip, 
the anisotropy term takes the form 

H A = -Nk z S 2 cos 2 ipnl - Nk z S 2 sin 2 ip n 2 x 

= -Nk z S 2 (cos 2 9 cos 2 ip + sin 2 6 sin 2 sin 2 ip) , (4) 

which describes simply the fact that components n z and 
n x profit from the easy axis as the ring is tilted around 
the field. 9 and <f> are the axial and azimuthal polar 
coordinates. The eigenvalues Ea are evaluated using the 
spherical harmonic wavefunctions Yi m (9, 4>). The results 
for the singlet and triplet levels are 

£o,o = -jA^S 2 (cos 2 ^ + sin 2 V) = ~Nk z S 2 , 

Ei,o = \Nk z S 2 cos 2 V - \Nk z S 2 sin 2 ^, (5) 

N 5 5 

E x ±x = — - -Nk z S 2 cos 2 ib - -Nk z S 2 sin 2 ib 
N 5 5 

±9VbB. 

Comparison of £0,0 with Ex,-x gives the result 



B c \ — 
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, ^iV 2 ^ 2 / 3. 2 , 

1 + — 1 - - s n 2 ib 

30J V 2 v 



(6) 



and thus a clear demonstration that the angle- 
dependence observed in both ferric wheel systems cor- 
responds to an easy-axis anisotropy for the spin direc- 
tion in the original model. For future reference, we 
define as B^ — AJ/g^sN the field value where the 
S = — > 1 transition would occur in the absence of 
magnetic anisotropy, whence Eq. (Q) may be reexpressed 
as 



B c i = B cl 



1 + — A c i 
15 



1 — — sin ip 



(7) 



where A c i is the value of the anisotropy-field raticcl 
A = 8k z JS 2 /(gfi B B) 2 at B = B%. Finally, by following 
the same procedure as in Eq. (|) for the wave function 
¥2,-2(9, <f>), one obtains 

12 7 1 3 

E 2 ,-2 = — - j Nk zS 2 cos 2 ip - -Nk z S 2 sin 2 ip 

-2gfi B B, (8) 
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and thus 



1+ 70 A 



(9) 



It is evident that in the presence of anisotropy, the widths 
of the magnetization steps are no longer equal. 

Before leaving the rigid-rotor model (Eq. (Q)), we may 
employ it to make a valuable qualitative observation con- 
cerning the nature of the magnetization. The situation 
of most interest for quantum transitions between low- 
energy states of similar energy is when the field is ap- 
plied perpendicular to the axis of the anisotropy. In 
the real materials, this anisotropy is rather weak, and 
it would be most appropriate to seek quantum coher- 
ent phenomena in the high-field regime of Ref. f|. If 
the field is dominant and sets the quantization axis 
[z] of the rotor, the on-site anisotropy term takes the 
form -J2i k z s x = -\Y Jl k z{S + + S-) 2 , where S a is 
a spin operator acting within the S — 5/2 manifold 
on each site. The anisotropy, considered as a pertur- 
bation, may change the total spin state of the system 
only by AS = 0, ±2. All level crossings between eigen- 
states of the rigid rotor will then be true crossings when 
Am = ±1, ±3, . . ., but will become anticrossings due to 
the mixing perturbation when Am = ±2, ±4, . . . Because 
the ground-state level alterations always involve Am = 1 , 
we thus obtain a simple understanding for the fact that 
all of the magnetization steps will have vertical portions 
for the starting Hamiltonian ([l]), whereas the energy level 
splitting Eoi — E\ — Eq wiU be not a sawtooth-shaped 
function but a rounded one.Q In principle, the anisotropy 
term may be reexpressed in terms of linear combinations 
of operators L a for the rigid rotor, allowing one to obtain 
analytical expressions for the numerical results to follow. 



III. EXPERIMENTAL COMPARISON 

The results of the previous section may be applied di- 
rectly to the experimental data of Refs. [| and ||. The 
variation of the critical fields for the lowest magnetization 
steps with the angle ip between the applied field and the 
ring axis has been measured to high accuracy by torque 
magnetometry, and is reproduced in Figs. 4 (Na:Fe@ and 
Feio) an d 3 (LiiFeg), respectively, of these works. The 
mean value of the field at which the step occurs gives 
the antiferromagnetic superexchange coupling J, which 
is responsible for maintaining the S = ground-state 
spin configuration at zero field. The variation of the crit- 
ical field values B c between their extrema at tp = and 
if) = 7r/2 gives the effective anisotropy parameter k z for 
each material. In the case of Li:Fe6, the two steps which 
have been characterized give an independent verification 
of the consistency of the rigid-rotor description over a 
wider field range. 
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Fe 10 


4.65T 


15.56K 


1.28T 


0.0088 
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Li:Fe6 


10.38T 


20.83K 


0.620T 


0.0053 


0.011 


Li:Fe6 


10.30T 


20.68K 


0.328T 


0.0066 


-0.091 



Fits of the data to Eqs. (||) or ([?]) for the first steps, 
and to Eq. (^|) for the second step in Li:Fe6, yield the 
values of J and k z shown in Table 1. The parameters 
in the last row are taken from B c i for Li:Feo, and show 
quite satisfactory consistency with the values deduced 
from the first step. The results are illustrated in Figs. 1 
and 2. The angular offset in the final column of the table 
is an alignment parameter concerning the experimental 
orientations of crystal and goniometer, and is not impor- 
tant for the intrinsic physics. This angle is quite large 
in Feio (Fig. 1(b)), where the crystal structure contains 
rings of two different orientations, which are observed^ 
in the same projection on the field direction. The fit- 
ting procedure employed here bypasses the definition of 
an intermediate, effective model with parameters Df, the 
results-.of Table 1 agree well with those of the original 
works,Ern which were extracted by different methods, and 
with the only other known torque measurements. □ They 
are also fully consistent with the energy scales deduced 
from analyses of low temperature specific heat for Na:Fe6 
and Feio,tJ and additionally for Li:Fe6.ll3 
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FIG. 1. Angular variation of critical field B c \ for first mag- 
netization step in (a) Na:Fe6 and (b) Feio. Data points are 
from Ref. and solid lines are fits using Eq. (Q) with the 
parameters in the Table. 
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where (S a ) — ^2iLi(Si. a ) is the ground-state expectation 
value of component a of the total-spin operator. The re- 
sults are shown in Fig. 3(a) for Fe6, with = 0.0136 J, 
and for angles ip — 0, 7r/4, and tt/2. The magnetization 
displays clearly the stepped form, with the step position 
varying as a function of field angle, which is observed in 
experiment and expected from the rigid-rotor model (|^) 
of Sec. II. From the location of the steps in field at each 
angle, we find that the model (Eqs. (0,^)) is in quantita- 
tive agreement with simulations at the 3-5% level. Given 
its inherent approximations, the rigid-rotor description 
is remarkably satisfactory in the experimentally relevant 
regime of weak anisotropy which is considered here. 
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FIG. 2. Angular variation of critical fields B c i (a) and B C 2 
(b) in Li:Fe6- Data points are from Ref. ^, and solid lines are 
fits using Eqs. (|^) and Q with the parameters in the Table. 



IV. NUMERICAL SIMULATIONS 

With this semi-classical background for the leading- 
order effects, we turn now to fully quantum mechanical 
simulations of the model (|l|) , with the minimal, uniaxial 
form of the U term 



N 



(10) 



We work with arbitrary field angle and strength, and with 
focus on the experimental parameter regimes for Na:Feg 
and Feio (Table 1). 
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FIG. 3. (a) Magnetization M(B) for Na:Fe 6 from ED, for 
field angles ip — 0, 7r/4, and tv/2. (b) Magnification of 
S — — > 1 step region, illustrating anisotropy-induced varia- 
tion in B c i, and vertical nature of step. 



A. Exact Diagonalization 

For Fe 6 , N = 6 spins S = 5/2 require dealing with 
(2S + 1) N = 46656 basis states. This Hamiltonian ma- 
trix is eminently accessible by the Lanczos exact diag- 
onalization (ED) technique, which provides the lowest 
eigenvalues and their corresponding eigenstates. For any 
angle of the field to the ring plane, the magnetization is 
given by 

M= — = {S z }cosiP + (S x )smil>, (11) 



The field resolution in these simulations was 0.001 at 
the steps in the units given. One observes that the mag- 
netization steps are vertical, as shown in Fig. 3(b) around 
the first step. However, for angles ip ^ they do not 
have the full height A M = 1. This feature, and the cur- 
vature of the magnetization "plateaus" , are illustrated 
more clearly in Fig. 4(a) for the same angles, but with 
the anisotropy taken as k z = 2fc^ a . They may be under- 
stood as a natural consequence of zero-point transitions 
between the nearly-degenerate energy levels of different 
spin sectors S z close to the step field B c \. These are 
enhanced over a broader region of field around any step 
B cn as the field angle approaches ip — tt/2, and the sys- 
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tern realizes more closely the optimal geometry.! of two 
degenerate spin states corresponding to the Neel vector 
n lying along ±£. Fig. 4(b) shows the magnetization 
curve for geometry ip = 7r/2, for three anisotropy values 
k z = 0, fc^ a , and 2fc^ a . It is clear that in the absence 
of the on-site anisotropy term no quantum transitions 
may occur, because all energy levels are degenerate with 
no barrier. Deviations of the critical fields from their 
unperturbed values B® n , and also the rounding of the 
plateaus, increase linearly with k z . Plateau raunding ef- 
fects are contained in the semiclassical theoryH, and may 
be compared at the qualitative level. 
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FIG. 4. Magnetization M{B) for an Fe6 ring from ED, (a) 
for field angles tjj = 0, 7r/4, and n/2 at anisotropy k z — 0.0272 
(twice that in Na:Fee), and (b) with field angle tp = tt/2 and 
anisotropy values k z = 0, fc^ a , and 2fc^ a . 



In Fig. 5(a) is shown the low-field torque , 

T y = -gfi B B((S x ) cos^j - (S z ) sm^j). 



(12) 



For the geometry ip = ir/A (illustrated), where the torque 
is maximal, it takes the form of a simple difference 
between the x and z spin components. As shown in 
Fig. 5(b), at higher fields this quantity is maximal around 
the 4th plateau, then decreases, changes sign, and adopts 
a negative constant value beyond saturation at the fif- 
teenth plateau. Such behavior may be understood from 
the asymmetric field evolution of the spin components 
in the presence of a finite easy-axis anisotropy benefiting 
one of them. The torque curve, which is directly experi- 
mentally accessible, is thus characterized by vertical steps 



at the magnetization jumps, doubly rounded plateaus, 
and a distinctive variation in step heights. The first two 
features may be expected to be visible at the lowest tem- 
peratures available in Fe6 systems, but the evolution of 
step heights would be hard to measure because of the 
high fields involved. 
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FIG. 5. Torque curve of Na:Fe6 for ip — n/4. (a) Detail at 
low fields; (b) full field range. 



B. Density Matrix Renormalization Group 

The DMRG technique^ is a powerful method for sim- 
ulating the low-energy properties of one-dimensional sys- 
tems, which can be much larger than those accessible by 
ED. Its extensions to periodic systems, and to high spins 
such as S — 5/2, are conceptually straightforward, but 
require considerable numerical effort due to the large spin 
degree of freedom and, for anisotropic ring systems in a 
field, the absence of symmetries. In our numerical calcu- 
lations, we retain at least 80 states, and use a finite-size 
algorithm with 3-5 sweeps for each case. With these con- 
ditions, the results are altered by less than 1 part in 10 3 
on adding more states or more sweeps. We present here 
the first results of simulations of this kind (Figs. 6-7). As 
for the ED results above, these are T = data, and show 
essentially the exact magnetization and torque response 
of 10-site rings for the minimal model of Eq. (|l|). Be- 
cause the DMRG method applied to a sufficiently small 
system (up to 8 sites) involves a very similar truncation 
to Lanczos ED, the results of simulations on Na:Fe6 are 
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identical to Figs. 3-5; this provides a useful check of both 
techniques, and obviates further discussion of this system 
by DMRG. 

Fig. 6 shows the magnetization curve of Feio for the 
physical parameter value k\° = 0.0088. The qualitative 
features of step location and plateau rounding with angle 
which are visible in the first three steps (Fig. 6(a)) are as 
in Figs. 3 and 4. Their relative sizes are similar to Na:Feg 
because the smaller value of k z is offset by the larger 
value of N. Indeed, the plateau rounding at ip = ir/2 is 
found to be more pronounced for Feio than for Na:Fe6, 
while, as emphasized above, the step remains vertical; 
the field increment in Fig. 6(a) is 0.0125 in the units 
displayed. Fig. 6(b) shows the magnetization response 
at high fields, which has 25 steps corresponding to the 
S z manifolds available between zero and full saturation. 
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FIG. 6. (a) Magnetization M(B) for Feio from DMRG, for 
field angles ip = 0, 7r/4, and ir/2. (b) M(B) for field angle 7r/4 
over the full field range, illustrating the 25 steps to saturation. 

Fig. 7 illustrates the torque for Feio at ^ = 7r/4. The 
qualitative similarity between Fig. 7(a) the data of Ref. ^ 
for the evolution of step heights is worthy of note, al- 
though the data at 0.45K have too much thermal broad- 
ening for meaningful comparison with the predictions for 
critical field locations and plateau rounding. However, 
because the critical field in this system is four times 
smaller than in Na:Fe6, the possibility arises of a de- 
tailed experimental characterization of this curve over 
the first five plateaus Following the torque response to 
high fields (Fig. 7(b)) shows that the maximal value ap- 



pears around the sixth plateau, and saturation after the 
25th. 



V. DISCUSSION 

The rigid-rotor model of Sec. II corresponds to the ki- 
netic limit, of a particle moving in a weak potential, and 
thus is quite suitable for the real Fe6 and Feio systems 
where the anisotropy k z is weak, and thus the barrier 
height is small. In addition to the energy levels and step 
positions discussed in Sec. Ill, the Hamiltonian (||) may 
be extended in second-order perturbation theory to give 
a qualitative description of the rounding of torque and 
magnetization plateaus. By contrast, the Lagrangean of 
Eq. (Q) is most appropriate for the tight-binding limit 
of a particle in a strong potential, and for the semi- 
classical description of tunneling dynamics. This formu- 
lation gives a useful qualitative account of the energy 
level splitting, and of tunneling processes which result 
in plateau rounding, but is not applicable to the experi- 
mental curves at low temperatures. The numerical sim- 
ulations remain the most reliable source of quantitative 
torque and magnetization data for the real system pa- 
rameters. 
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FIG. 7. Torque curve of Feio for ip = 7r/4. (a) Detail at 
low fields; (b) full field range. 

With regard to the possible observation of quantum 
coherent tunneling of the Neel vector in ferric wheel sys- 
tems, we may state that the magnetization and torque 
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curves do appear-, to demonstrate that the desirable 
physical situationu is realized. This consists of two, 
nearly-degenerate, low-lying energy levels, well separated 
from higher levels and weakly mixed by quantum (zero- 
point) transitions through a barrier, due to the easy- 
axis anisotropy, which separates them. However, the key 
property which cannot be shown by thermodynamic mea- 
surements such as the magnetization is that of coherence, 
or evolution of the quantum state over several oscillation 
cycles without scattering into a different state (decoher- 
ence). The experimental and theoretical aspects of this 
question, and the quantities (time-dependent correlation 
functions) which may be computed analytically and nu- 
merically to provide additional insight, will be addressed 
in a future publication.E3 

We conclude our analysis with a brief discussion of its 
relation to experiment. The existence of magnetization 
and torque steps, their locations for different field ge- 
ometries, and their heights, are in good agreement with 
torque magnetometry results obtainedtj for Fe6 and Feio 
samples at temperatures of 0.45K. We stress again that 
the simulations for the minimal model of Eq. ( |l0| ) show 
a vertical step at the appropriate critical field for all 
field angks. However, the recent, highly accurate results 
obtainedH for the first torque step in Na:Fe6 at temper- 
atures down to 30mK show two important features not 
reproduced by simulations of the model. One is the lin- 
ear torque response observed for the 5 = plateau, and 
the other a temperature-independent field broadening of 
the magnetization steps. 

The first feature may be caused by admixtures of 
higher-spin states in the composition of the ground state, 
of either "extrinsic" , meaning crystal misalignment or 
twinning, or "intrinsic" origin, meaning additional inter- 
action terms not contained in (|l0|). Of yet more interest 
is the second feature. An intrinsic origin for this effect 
implies a term in the Hamiltonian permitting transitions 
of AS Z = ±1. Such terms do not arise from biaxial 
anisotropy, particularly in a ring geometry. This in turn 
may indicate a possible role for superhyperfine interac- 
tions with proton spins on H-sites neighboring the metal- 
lic ions. However, a studyli3 of this interaction in Feio 
has found that its effects are too small by some orders 
of magnitude to explain the observed broadening. Sim- 
ilarly, hyperfine interactions with the Fe nuclear spins 
may be excluded because the natural abundance of 57 Fe, 
the only species with non-zero nuclear spin, is so low. 
The most realistic intrinsic possibility is that the step 
width represents the energy scale of inter-molecular inter- 
actions, as any coupling between the moments of neigh- 
boring rings would broaden the magnetization response. 
The observed full width of the first step at 40mKH is 
0.4T, which would correspond to an effective interaction 
of 0.14K per Fe atom in each ring. While the possible 
contributions to such interactions are difficult to quantify, 
simple geometricaLeonsiderations based on the crystal- 
lographic structures suggest that the long-range dipolar 
component could indeed be of this order. 



A candidate extrinsic originH is that molecules at the 
crystal surface may have different lattice parameters due 
to solvent loss, and thus undergo transitions at different 
field values corresponding to altered interaction strengths 
J. This effect should lead to rather small tails in the 
distribution of critical fields compared to the bulk re- 
sponse. Crystal mosaicity presents a further extrinsic 
possibility for step broadening, in that molecules with 
different angles to the applied field undergo transitions 
at different critical fields. In this scenario, the full width 
of the broadened first step at 40mK would correspond to 
a quasi-Gaussian angular distribution with a full width 
of 9.2°, a rather large number which should be readily 
falsifiable from x-ray structural data. A more detailed 
analysis of these possibilities is deferred until further ex- 
perimental data emerge to clarify the issue of extrinsic 
or intrinsic origin. 

VI. CONCLUSION 

We have presented a detailed analysis of the magneti- 
zation and torque response of the molecular magnetic sys- 
tems Feg and Feio- The step structure with critical fields 
B cn , and variation of B cn with field angle to the ring 
plane, are readily understood from low-spin subspaces 
of a rigid-rotor Hamiltonian with weak spin anisotropy. 
Numerical simulations of the minimal but fully quan- 
tum mechanical model demonstrate clearly a vanishing 
step width, accompanied by rounding of the plateau re- 
gions even at zero temperature. These features can be 
related to those zero-point transitions permitted within 
this model when the system is close to magnetization 
steps, where energy levels of different spin states are in 
close proximity. The differences .between simulation and 
experimental observation in FegB demonstrate the need 
for additional physical considerations on the theoretical 
side, and characterization of a wider range of samples on 
the experimental, with a view to discerning those features 
truly intrinsic to the ferric wheel systems. The angular 
variation of the step and rounding effects, which emerge 
from simulations and theory, should be observable in 
the most accurate experiments currently possible.El These 
would be of particular interest in Feio, and may be ex- 
pected to verify that the ferric wheel systems in trans- 
verse magnetic fields represent good candidates for the 
observation of quantum coherent tunneling of the Neel 
vector. 
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